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Abstract 

Two microscopic theories for multiple resonances in nuclei are compared, 
n-particle-hole RPA and quantized Time-Dependent Hartree-Fock (TDHF). 

The Lipkin-Meshkov-Glick model is used as test case. We find that quantized 
TDHF is superior in many respects, except for very small systems. 

1 Introduction 

The most prominent excitations of the nucleus are the giant resonances. One has 
therefore dreamed about multiple excitation of giant resonances quite early in the 
development of nuclear theory [|I|. Recent experimental progress has brought those 
multiple resonances into the range of observability [^. This has inspired several 
theoretical investigations in the framework of a second RPA, see e.g. |P, 
Multiple excitations of basic collective modes are well known from the low energy- 
spectra of soft nuclei. They are most often characterized phenomenologically in 
terms of interacting elementary Bosons, e.g. surface vibrations and pairing modes 
PJ. Theoretical foundations for such a collective description of the low-energy modes 
have been long searched for in the framework of the Boson expansions, for a review 
see [0. This task has turned out to be very involved, and a possible explanation is 
that the low energy modes are far from a simple picture of collective motion as mere 


1 


shape vibrations of the mean held. Motion along a deformation path is accompanied 
by dramatic reoccupations amongst crossing single particle levels [|] having possibly 
curious consequences for the collective dynamics 0. Giant resonances, on the other 
hand, represent the much cleaner collective modes, dominated by the mean held and 
without noticeable pairing ehects. The theoretical appeal of the multiple resonances 
is then to study large amplitude motion in a case where the underlying elementary 
Boson is fairly well understood. The basic modes, the giant resonances, are usually 
described within the RPA, the theory of small Iph excitation about the Hartree-Fock 
ground state []^. Multiple resonances require extensions beyond standard RPA. 
A straightforward next step is the ’’second RPA” extending the treatment to the 
space of 2-ph excitations Boson expansions, which have also been considered 

in that context conhne the extension to the collective modes but allow to go 
further in the anharmonic terms. Somewhat more general are the ’’higher RPA” 
approaches considering coherent mixtures of n-ph-excitations. The equations-of- 
motion techniques as outlined in 0,0 provide an extremely useful tool to derive 
those expansions along a given set of n-ph operators. A much different point of view 
is presented by Time-Dependent Hartree-Fock (TDHF) approach which pronounces 
more the vibrating mean field. Its small amplitude limit delivers also the RPA. 
The spectra for higher modes can be computed with quantized TDHF which picks 
the quantum states amongst all the classically allowed energies by requiring strictly 
periodic TDHF orbit which have integer action along one closed orbit El- It is 
the aim of this paper to compare the results on multiple resonances as computed 
with quantized TDHF or with an n-ph-KPA derived from an equations of motion 
technique. The comparison is performed in a schematic model for nuclear dynamics 
in the active shell around the Fermi surface, the Lipkin-Meshkov-Glick (LMG) model 
BH- We will use in the following the abbreviation RPA-At for such a higher RPA 
which includes up to N-ph excitations. 

The paper is outlined as follows: In section |^, we introduce briefly the LMG 
model. In section we recapitulate the equations-of-motion technique and work 
it out for the present model. In section we explain quantized TDHF and how 
we compute it. Finally in section |^, we present and discuss the results from both 


schemes. 


2 




2 The Lipkin Model 


The LMG model [[l^ |T^ is a simplified shell model with two degenerate bands of 
single Fermion states s G {+,—}, m G The uncoupled ground state 

of the model carries all states s = — as occupied and s = + as empty. Due to 
the degeneracy, all single particle excitation energies are the same e. The model 
Hamiltonian reads 


H = eK, + x ‘-^A K^K+ + K.K.) (1) 

where the basic operators are 


Ko 

— r 1 1 





K+ 


(2) 




K_ 

= A 

(3) 


The second term, quadratic in the K±, models the residual interaction. The K± are 
Iph operators with respect the uncoupled ground state and thus the Hamiltonian 
(|^) is tailored for model studies of the dynamics along the Iph channel, i.e. for mean 
field dynamics. The dimensionless parameter y describes the interaction strength 
to the single=particle splitting e. The scaling with 1/(^2 — 1) places the critical 
coupling at X = 1 almost independent of the system size D. 

Note the two important features of the LMG model; First, the residual interac¬ 
tion is only active for vertical Iph excitations, i.e. those which leave the secondary 
quantum number m invariant. Second, all substates m are handled with the same 
weight and phase in each of the three operators Ka- This introduces a high symme¬ 
try in the model such that we have an exactly decoupling collective subspace which 
is spanned by the Thouless transformations |<F) = exp (^2 X]Qe{-,o,-i-} |0) with 

the basic collective operators . Moreover, these basis operators form a simple 

SU{2) algebra 

[Ko,K±] = ±K± and [K+,K_]=2Ko . (4) 

All further evaluations, in the equations-of-motion technique as well as in quantized 
TDHF, require only this algebra and, of course, the model Hamiltonian (|^. 


3 










In spite of its simplicity, the LMG model unfolds a rich variety of scenarios for the 
collective motion, nearly harmonic motion for small Xi increasing anharmonicities 
with increasing y, and unstable motion with subsequent phase transition near y = 1, 
for a visualization see the beginning of section It provides the ideal testing ground 
for our purposes. In fact, it was one of the originals aims of the LMG model to 
understand the appearance of collective modes in the RPA. It is the mode which is 
excited by a mix xKj^ — yK^ and whose frequency is a; = ca/I ~ in the domain 
X < 1- 

3 Equations-of-motion and RPA-A^ 

The equations-of-motion technique as outlined in Jll], ^ is in our opinion the most 
obvious way to formulate a RPA-iV as an algebra of selected excitation operators. 
The aim is to optimize a set of excitations operators Qj, which excite the state 

w = oiio) (5) 

out of a correlated ground state |0). The excitation operators are taken as linear 
superposition from a basis set of operators Qa 

= . ( 6 ) 

OL 

In the ideal case, the excitation operators should fulhll he Heisenberg equation 

[H,Ql] = {E,-Eo)Ql . (7) 

But this can rarely be fulhlled in practice because the algebra of the and El is not 
closed. One thus requires Eq. (|^) in the average after ’’stabilizing” the expression 
to a double commutator. This leads to the variational equations 

<0\[SQ,[H,Ql]]\0>=uj,<0\[6Q,Ql]\iy> , oo, = {E, - Eo) ( 8 ) 

where 6Q stands for any variation in the given space of Qa- The correlated ground 
state |0) can be generated with a similar variational equation. It is written in the 
form 

|0>=e^|i7F> (9) 
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where S is an antihermitian operator and \HF > the Hartree-Fock gronnd state. 


The Operator S is determined by [11 


< 0\[6S, i7]|0 >= 0 


( 10 ) 


where SS is again any variation in the given basis space of operators. 

The spectrum of collective excitations in the Lipkin Model is generated from the 
basis states iF? and K^. We thus make for the excitation operator the ansatz 


N 


Ql = J2xn,uKl-yn,uK^ 


( 11 ) 


n=l 


and similarly for the correlation operator 

M 


s = Y.^n,u {ki-k: 


n=l 


( 12 ) 


The N and M represent the order of the expansion, in practice the maximal order 
of N-ph or M-ph operator taken into account. The equations-of-motion and 
constitute a set of equations to determine cUy, Xi^u, yi,u {n = l..A^, v = l..iV) 
and 7 „ (n = 1..M). It is most consistent to use the same order M k. N in both 
pieces and this yields usually the best results for a given expense. An exception is 
the conventional RPA which appears as the case N = 1 and M = 0 within that 
scheme. The enhancement to M = 1 is unnecessary in that case because the Iph 
states decouple from the Hartree-Fock ground state []^. For any other A^ > 1, it is 
best to chose N = M and this is what we call the higher RPA of order N^ in short 
RPA-A^. We will consider in the following the two particular examples RPA-2 and 
RPA-5. 


4 Quantized TDHF 

TDHF approximates the dynamics of a many-Fermion system by a time-dependent 
Slater state |*h(t)). The time-evolution of |*F(t)) is optimized by deriving the TDHF- 
equations from the variational principle of stationary action 

5 ( <^ih^-H\^> dt = Q (13) 

J at 
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where 6 denotes variations within the snbspace of Slater states. TDHF delivers a 
deterministic eqnation-of-motion for |*h(t)). It can be interpreted as representing 
the classical limit of many-Fermion dynamics. The interpretation is pertinent in 
particnlar with respect to the fact that TDHF motion is possible for any initial 
condition |*h(0)) at any energy. However, it can be complemented by a (semiclassical) 
qnantization, as has been shown by very different approaches, as e.g. the fnnctional 
integral representation of the Green fnnction G{E) = tr--^ |]^ |^, 


or the 


method of gange invariant periodic quantization The quantum states are 

those solutions of the TDHF equation (0 which fulhll two additional conditions: 
First, they are strictly periodic including the phase. 


|4(( + T)> = |4(()) , 

and second, they fulhll the quantization condition 

/ < <hhh—|<F >= 2n7rh 

Jo dt' 


(14) 


(15) 


where n is some integer number. 

The most problematic part of quantized TDHF is to hnd periodic solutions 


in the chaotic manifold of trajectories of a realistic TDHF calculation. Fortunately, 
this is no problem at all in the LMG model because the motion of interest is reduced 
to the one collective degree-of-freedom spanned by the K± as dehned in section ||. 
First tests of quantized TDHF have therefore been performed within this model 
0, 1^ . According to the Thouless theorem, each Slater determinant of a given 


system can be generated from a reference determinant |<ho) as |$) = exp (Ap/i)|<Fo) 
where Aph is a Iph operator. We aim to span the space of ” collective deformations” 
outgoing from the Hartree-Fock ground state. This is achieved by the coherent 
states 


\z >= e 


zK^ 


\HF> {l + \z\Y^ 


(16) 


which are labelled by the complex shift parameter z. Note that these coherent states 
are already normalized. 


The expectation values which are necessary for Eq. ([1^) are easily evaluated 
using the algebra (§) and the subsequent simple properties of the coherent states. 
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We obtain 


(17) 

(18) 


< z\i%—\z > 

' dr 

< z\H\z > 
z 


— — h{l — cos-0)0 
— Y (cos '0 ~ ^ sin^ -i/i cos 2(j)) 


(19) 


The time evolution of the coherent state is contained in the time-evolution of the pa¬ 
rameter By Variation according to (P!B|) we hnd the TDHF equations as equations- 
of-motion for z 



e 



e 


—X sin d sin 20 
— (x cos 0 cos 20 -|- 1) 


( 20 ) 

( 21 ) 


The periodicity condition (JT^) is trivially fulhlled because we have an effectively 
two-dimensional phase space. It remains to watch the quantization condition 
which reads in the LMG model 


n 


(■T . 

(1 — COS 0)0 df =- 


1 — COS0) dcj) = 2n7r 


The initial condition of the coupled differential equations 


and 


( 22 ) 


are varied 


until the solution fulhlls the quantization condition (|^) . The (quantum) excitation 


energy is then provided by the expectation value ([Tq) which is constant along the 
stationary path. 


5 Results 

The LMG model has essentially two parameters, namely X; the relative coupling 
strength, and G, the particle number or size of active phase space. Variation of x 
varies the ground state deformation and the anharmonicity of the collective reso¬ 
nances. The uncoupled ground state of the model is the Hartree-Fock ground state 
for X < 1 whereas a transition to a ’’deformed” state occurs for x > 1- Glose to har¬ 
monic excitations appear at x 1 and x ^ 1 whereas increasing anharmonicities 
build up if the critical point x = 1 is approached from both sides. Variation of G 
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influences the softness of the transition around the critical x = being softer for 
small n and developing a sharp phase transition for hi —> oo. In order to visual¬ 
ize the dynamics of the system, we show in Fig. |I] the deformation energies versus 
collective deformation angle V’- 

The case y = 0.4 displays a nearly parabolic potential indicating that we have 
to expect essentially harmonic motion there with small anharmonic perturbations. 
The case y = 1 shows an extremely soft and anharmonic potential. The system 
is in a critical regime where small perturbations lead to huge reactions. The ef¬ 
fect remains finite for finite hi and tends to critical, i.e. infinite, fluctuations for 
n —>• cx). The case y = 4 shows again a well developed ground state with small 
oscillations about it in the first excited state. But note that this ground state is 
placed at finite deformation and that it appears twice at exactly the same energy. 
As a consequence, we will have a doubled excitation spectrum because each mode 
has a symmetric copy about the other ground state. The excitation energies from 
quantized TDHF or RPA-A^ are exactly degenerate. The degeneracy of the exact 
excited states is removed a little bit due to tunneling processes through the barrier 
at -0 = 0. None of the presently discussed theories can reproduce this splitting. 

In Fig. ^ we compare the results of quantized TDHF, RPA-5 and RPA-2 with 
the exact solution up the fourth resonance excitation. The figure concentrates on 
the regime of undercritical couplings. Three system sizes are considered which are 
comparable in the nuclear case to the phase spaces of light, medium and heavy 
nuclei. The first excited state is very well described by all methods. The tiny 
differences seen there would favour quantized TDHF and, somewhat surprisingly, 
RPA-5 looks a bit inferior to RPA-2 for small D. But one should not overstress 
these details. The comforting message is that plain RPA performs well for the first 
excited state in the regime of undercritical coupling. More visible differences show 
up with increasing excitations. Quantized TDHF looks inferior for the small system 
size D = 5 where RPA-5 performs very well. The reverse happens for medium and 
large systems where the result from quantized TDHF are astonishingly close to the 
exact solution. This complies with the interpretation of TDHF as a classical limit of 
the many-Fermion dynamics. This is a classical limit in the sense of a 1/D expansion 
and it becomes increasingly valid with increasing system size, here D. On the other 
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hand, the success of the RPA-5 for f2 = 5 is not so surprising because that expansion 
is then close to complete (note that it is not yet fully complete at 5th order because 
the 0-ph operator Kq is missing in the expansion and has to be regenerated from 
pairs of K±). Finally, we note that the RPA-2 can by dehnition only describe the 
system up to the second excited state. But it performs surprisingly well at its upper 
limit UJ 2 and can reliably be used instead of RPA-5 there. 

In Fig. ^ we show the same comparison for a wider range of couplings y. The 
exact spectra reflect nicely the appearance of two equivalent ground states in that 
pairs of excited states develop for y — oo. The hrst excited state at y < 1 goes to 
zero excitation and thus merges into one of the two ground states of the system, in 
fact it represents the antisymmetric combination whereas the former ground state 
develops into the symmetric combination. The transition is best developed for the 
large whereas the small = 5 stays far from the asymptotic stage within the 
range of y considered. It is interesting to see that the different approximations 
handle the transitional region much differently. The plain RPA is known to break 


down near the critical y This feature persists also for the higher RPA-A^. 

Even the RPA-5 runs into difficulties if y grows near 1. This region has to be 
excluded from any RPA-like expansion. It is a typical application case for theories 
of large amplitude collective motion, as e.g. the generator-coordinate method 


RPA becomes applicable again beyond the transition where well bound ground state 
minima develop again, see Fig. |l|. Accordingly, a new branch of RPA solutions can 
be computed. We show in Fig. the results from the RPA-2 because we have there at 
most a second excited state and we know from Fig. ^that this level of approximation 
is fully sufficient up the second excitation. The RPA solutions in the regime of large 
y look fair for large y but run also in insurmountable difficulties if y is lowered 
towards the critical y = 1. We thus have a separate branch of RPA solutions for 
large y which is disconnected from the branch of RPA solutions in the regime of 
small couplings. As expected from Fig. |l|, the motion in the deformed minima is 
much less harmonic. As a consequence, the RPA-2 gives less reliable results in the 
regime y > 1, in particular for the small system. For example, it misses the coupling 
through the barrier and gives always exactly degenerate pairs of excitation states. 

It is interesting to see that quantized TDHF is much more robust at the critical 
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point and provides surprisingly appropriate results there. It is indeed able to stretch 
branches of solutions over all ranges of the coupling y. However, a bit of care is still 
required because once in a while a branch changes discontinously. This happens 
in the strong coupling regime whenever the quantized energy crosses the barrier 
between the left and the right minimum. TDHF as a classical theory cannot describe 
tunneling through the barrier. Thus the branch of extended motion over the barrier 
ceases to exist and two new branches below the barrier develop in the separate 
minima. To give an example, the first excited TDHF state in case of D = 5 stops 
at X ~ 3. This reflects the fact that it should merge into the second ground state of 
the system, see the exact solution for comparison. Thus the dashed line in the hgure 
should thought to be discontinously connected to ci;i = 0 and this behaviour serves as 
an approximation to the exact curve, a crude approximation indeed, but at least the 
behaviour is qualitatively present. Similar considerations apply in all the other cases 
where a TDHF branch stops. We thus see that quantized TDHF gives a very robust 
description in all ranges of y. It is quantitatively correct almost everywhere, except 
at the discontinuities of the branches. It provides at least the proper qualitative 
behaviour in those critical cases where branches of solutions disappear or creep up. 
And, these critical cases can be clearly identihed such that one knows where to trust 
the solution and where better not. 


6 Conclusion 

In this paper, we have investigated two microscopic theories for multiple resonance 
excitations in nuclei, the RPA-iV, an extension of the RPA to higher n-ph states, 
and quantized TDHF. The RPA-iV stands for the class of theories which are based 
on opertor expanions and the Heisenberg equations-of-motion. A close relative to 
RPA-iV are the Boson Expansions such that our Endings apply for these as well. 
The RPA-iV is an expansion in orders of excitations which is easiest and most reli¬ 
able for low excitation modes and becomes increasingly cumbersome with increasing 
energies and increasing particle number. The quantized TDHF, on the other hand, 
concentrates on the time-dependent state of the system in the Schrbdinger picture. 
The mean-field approximation, employed in TDHF, corresponds to the classical limit 
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of the dynamics of the many-body system. It thus becomes increasingly valid with 
increasing energy and increasing particle number. From this formal point-of-view, 
RPA-A^ and quantized TDHF cover contrary regimes. RPA-iV starts an expansion 
from the lowest excitation modes and stays fully quantum mechanical throughout 
each step whereas quantized TDHF as a classical approach looks more appropriate 
for high energies and has to be semiclassically requantized. 

The practical comparison has been performed within the Lipkin-Meshkov-Glick 
(LMG) model which mimics the nuclear excitation dynamics as well as a structural 
transition from symmetric to deformed ground states. We have studied three typical 
cases representing a small, a medium, and a heavy nucleus. The RPA-iV is found 
to work hne in the domain of well bound systems, i.e. far below and far above the 
critical transition point. However, it runs into serious difficulties around the critical 
point which cannot be surmounted even by high order expansion (actually we have 
tested it up to N=5). Quantized TDHF, on the other hand, provides surprisingly 
robust results for each excitation energy and in each structural domain including 
the transition point. The classical character of TDHF shines through in the fact 
that it is somewhat less reliable for the smallest system in our sample. We thus 
conclude that quantized TDHF is superior in most domains, but RPA-iV remains 
the method of choice for small systems at not too high excitation energies. 

The simplicity of the LMG model inhibits practical estimates of the expense 
of the methods. The RPA-iV at full extend will become unfeasible for a realistic 
nucleus, except perhaps N = 2 for small nuclei. Higher N can be treated as Boson 
Expansion with a basis of fixed 1-ph excitation operators. The necessary algebra of 
multiple commutators can be carried through nowadays but will become very cum¬ 
bersome. Quantized TDHF, on the other hand, is formally very simple as it employs 
always as a unique ingredient a TDHF code which is nowadys easily available. But 
it poses the serious technical problem to find the branches of periodic orbits in the 
chaotic multitude of solutions. This could turn out to be a big hindrance for the 
application of the otherwise very straightforward and robust approach. 

We would like to thank G. Toepffer for his constant interest on the work and many 
helpful comments. D.K. and M.W. have been supported by the Studienstiftung des 
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Figure 1: 

Visualization of the deformation energy (solid lines) versus the collective deforma¬ 
tion angle ip for the LMG model with = 14 particles at three different coupling 
strengths as indicated. The energy of the hrst excited mode is indicated by a dashed 
line in each case. 


Figure 2: 

Spectra of the LMG model as a function of the strength of the residual interaction in 
the range of undercritical coupling x < 1 for three particle numbers as indicated. 
The different line types represent different approximation schemes: solid = exact 
solution, dashed = quantized TDHF, dash-dotted = RPA-5, dotted = RPA-2. 


Figure 3: 

Spectra of the LMG model as a function of the strength of the residual interaction 
in a broad range of couplings for three particles number Gas indicated. The different 
line types represent different approximation schemes: solid = exact solution, dashed 
= quantized TDHF, dash-dotted = RPA-5, ’spherical’ ground state and RPA-2, 
’deformed’ ground states, dotted = RPA-2, ’spherical’ ground state. Note that a 
RPA-5 on a ’spherical’ ground state and a RPA-2, ’deformed’ ground states both 
describe the hrst six eigenstates of the Hamiltonian. 
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